A distributed architecture for scalable quantum computation with realistically noisy 

devices 
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Tremendous efforts tiave been paid for realization of fault-tolerant quantum computation so far. 
However, preexisting fault-tolerant schemes assume that a lot of qubits live together in a single 
quantum system, which is incompatible with actual situations of experiment. Here we propose a 
novel architecture for practically scalable quantum computation, where quantum computation is 
distributed over small-size (four-qubit) local systems, which are connected by quantum channels. 
We show that the proposed architecture works even with the error probability 0.1% of local oper- 
ations, which breaks through the consensus by one order of magnitude. Furthermore, the fidelity 
of quantum channels can be very low ^ 0.7, which substantially relaxes the difficulty of scaling- 
up the architecture. AH key elements and their accuracy required for the present architecture are 
within reach of current technology. The present architecture allows us to achieve efficient scaling of 
quantum computer, as has been achieved in today's classical computer. 



I. INTRODUCTION 

Scalability and fault-tolerance are essential ingredients 
in both classical and quantum computation. In 1945, 
von Neumann proposed the first architecture design of 
a classical computer and most of today's classical 
computers are of von Neumann type. He also explored 
fault-tolerance in classical computation by developing a 
systematic way to construct reliable logic circuits from 
noisy devices 2* , which is now attracting renewed interest 
in the field of nanocomputer Q . 

In quantum computation, on the other hand, it is still 
under extensive investigation what type of physical sys- 
tem is best suited to an experimental realization of quan- 
tum computation and what level of accuracy is required 
for it. For a better understanding of them, a lot of efforts 
have been made over the past decade both theoretically 
and experimentally. Quantum fault-tolerance theory en- 
sures scalable quantum computation with noisy quantum 
devices as long as the error probability of such devices 
is smaller than a threshold value (see Ref. [1] and refer- 
ences therein) . The noise thresholds have been calculated 
to be about ^ 0.1-1% for several fault-tolerant schemes 
under various assumptions f^-Q- In the experiments, it 
is nowadays possible to control a few to dozen qubits in a 
wide variety of physical systems such as trapped ions [1| 
and nitrogen-vacancy (NV) centres in diamond @ etc. 
(see Ref. ^Oj and references therein). There is, however, 
still a large gap between these top-down (theoretical) and 
bottom-up (experimental) approaches. For example, the 
existing fault-tolerant schemes [5-7] require many qubits 
live together in a single system and assume that the whole 
system can be controlled with the same accuracy regard- 
less of its size. In experiment, on the other hand, the 
number of qubits in a single system is rather limited; if 
we increase the number of qubits in a single system, the 
control becomes more and more complex, which makes 
it hard to achieve the same accuracy. In fact, as men- 
tioned in Ref. ^lli] there is a consensus that "for practical 



scalability the probability of error introduced by the ap- 
plication of a quantum gate must be less than 0.0001". 
In order to fill the gap and break through the consensus, 
we have to develop a fundamental design of a scalable 
and fault-tolerant architecture for quantum computation 
as von Neumann did in the early years of classical com- 
puter technology. 

Several architectures for scalable quantum computa- 
tion have been proposed so far in bottom-up approaches 
by considering specific physical implementations such as 
photons [12] and semiconductor nanophotonics [l^. It 
is, however, still not fully understood what kind of ele- 
ments and what level of accuracies of them are essentially 
required to construct a scalable and fault-tolerant archi- 
tecture. Furthermore, the requirements of the existing 
quantum architectures are too demanding and compli- 
cated to be experimentally feasible. In this paper, we 
propose a fundamental and simple design of a scalable 
and fault-tolerant architecture for quantum computation 
taking a top-down approach, and clarify under what con- 
dition scalable quantum computation can be executed 
fault-tolerantly. 

To ensure scalability, we adopt a distributed approach 
to quantum computation [l4j , where small quantum sys- 
tems, say quantum arithmetic logic units (QALUs), are 
connected via quantum channels. This type of archi- 
tectures are equipped with built-in modular scalability. 
That is, the size of the architecture scales up by adding 
well-established individual units. In the pre vious one- 
dimensional (ID) distributed architecture [Tsl. IT6|. how- 
ever, entangled states have to be shared between QALUs 
of arbitrary distance, and therefore the quantum repeater 
protocol [14] is utilised. As the distance gets larger, this 
procedure takes more time. Moreover, every QALU in- 
between must provide a workspace for this procedure, 
resulting in an overhead which grows with the distance. 
Besides, if there is a single point of breakdown in the 
quantum channels possibly due to a manufacturing er- 
ror, one cannot generate any entanglement between two 
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FIG. 1: The 2D architecture and QALUs. The 2D dis- 
tributed architecture, which consists of four-qubit (blue cir- 
cles) QALUs (boxes) connected by quantum channels (lines) 
with their nearest-neighbours. 



parts separated by this point. 

The present architecture, on the other hand, consists 
of a two-dimensional (2D) array of QALUs, each of which 
consists of four qubits, as depicted in Fig. [1] Fur- 
thermore, entangled states are shared only between the 
nearest-neighbour QALUs, which allows highly parallel 
operations. Surprisingly, this simple architecture is found 
to be sufficient for fault-tolerant quantum computation 
and works with very noisy quantum channels of fidelity 
~ 0.7 and reasonably accurate local operations of the 
error probability ^ 0.1%, which breaks through the con- 
sensus by one order of magnitude These results 
are achieved by utilizing twofold error management tech- 
niques: entanglement purification |l8l42C| and topologi- 
cal quantum computation (TQC) [6[. The former is em- 
ployed to implement a reliable two-qubit gate by using 
very noisy quantum channels with the help of quantum 
gate teleportation In particular, we apply high- 
performance entanglement purification, so-called double 
selection scheme |2C| , which is essential for achieving the 
above result. TQC is used to handle the remaining er- 
rors and to archive quantum gate operations of arbi- 
trary accuracy, which is required for large-scale quantum 
computation. Furthermore, the nearest-neighbour quan- 
tum communications together with TQC is quite robust 
against the manufacturing error mentioned above, since 
we can reconstruct a reliable logical information on the 
surface code by avoiding such defects [H, [l^j • 



All key ingredients in the present architecture, (i) a 
four-qubit system, (ii) gate operations in the four-qubit 
system, and (iii) entangling operations between the sep- 
arate systems, have already been demonstrated exper- 
imentally in various physical systems. Actually, the 
benchmarks in trapped ion systems are comparable to the 
requirements of the proposed architecture. These results 
push the realization of large-scale quantum computation 
within reach of current technology. We believe that this 
work providing a fundamental and simple design of a scal- 
able and fault-tolerant architecture for quantum compu- 
tation fills the gap between the top-down and bottom-up 
approaches, and gives a good guideline and benchmark 
in development of devices for quantum computer. 



II. 2D DISTRIBUTED ARCHITECTURE 

The present architecture consists of a 2D array of 
well-defined small (four-qubit) local quantum systems, 
QALUs (see Fig. [T]), where we can implement the ini- 
tializations, measurements, and two-qubit gate opera- 
tions. These local operations are assumed to be im- 
perfect, which is modeled as follows: (i) An ideal two- 
qubit gate is followed by two-qubit depolarizing noise, 
(l~Ps)P + E(j,i)^(o,o)Kj(o-*<^o-j)p(o-i(8)crj) where Pg = 
X](i,j)#(o,o)Py = 0,1,2,3), and tr^ {i = 0,1,2,3) are 
the Pauli matrices (ctq = I)- (ii) The measurement of a 
physical qubit is implemented with an error probability 
Pm- The QALUs are assumed to work with reasonably 
high accuracy {pg,PM ~ 0.1%). The memory error prob- 
ability is assumed to be sufficiently smaller than that of 
the operational errors (as is also the case for most of phys- 
ical systems) and ignored for clarity (otherwise, since the 
number of waiting steps I is finite in the present 2D archi- 
tecture, we can take the memory errors into account by 
replacing the error probability Pg of the two-qubit gates 
with Pg + r]l, where i] indicates the memory error proba- 
bility per step) . Such QALUs are connected via quantum 
channels of fidelity F. Here the channel fidelity F means 
that we can share a maximally entangled state (MES) of 
fidelity F between the nearest-neighbour QALUs. The 
remote entangling operations can be used for this pur- 
pose j23| . The quantum channels are assumed to be rela- 
tively noisy {F ^ 0.7-0.9), which relaxes the complexity 
in scaling up the size of the architecture. 



III. ENTANGLEMENT PUMPING AND 
TELEPORTATION-BASED TWO-QUBIT GATE 

The noise introduced by the quantum channel degrades 
the accuracy of the gate operation between the separate 
systems. In classical computation, this kind of obsta- 
cles can be overcome by using the multiplexing technique 
where multiple copies of data with majority voting 
are used to ensure the reliability of logic circuits. In the 
case of quantum physics, however, we cannot make copies 
of an unknown quantum state due to the no-cloning the- 
Instead, we can use entanglement purifica- 
igj and quantum teleportation [2l| in order to 
improve the fidelity of the gate operations between two 
QALUs. In particular, a novel method is devised to re- 
alise gate operations of high accuracy, where an intelli- 
gent use of the four-qubit local system allows us to obtain 
MES of high fidelity even with very noisy quantum chan- 
nels, as shown below. 

The noisy copies of the MES (we call it "noisy MES" 
hereafter) , which are shared by using the noisy quantum 
channels, are purified by using the entanglement pump- 
ing scheme [lajISal as follows [see Fig. [11(c)]: (i) A noisy 
MES (we call it the target pair) is purified by reducing its 
bit-flip (A) error via level- 1 entanglement pumping with 
double selection [23|, where two noisy MESs are used as 
the ancillae. The phase-flip (Z) error in the target pair 
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FIG. 2: Entanglement purification with entanglement pump- 
ing, (a ) The entanglement purification with single selection 
[la . . (b) The entanglement purification with double se- 
lection [20]. (c) (i) The level-1 entanglement pumping with 
double selection, where one noisy MES is purified by using 
two noisy MESs as ancillae. (ii) The level-1 entanglement 
pumping with single selection, where one noisy MES is puri- 
fied by using one noisy MES. (iii) The level-2 entanglement 
pumping with double selection, where the output state of (i) 
is purified by using the output state of (ii) and one noisy MES 
as ancilla. 



increases in this process, (ii) Another noisy MES is pu- 
rified by usin g le vel-1 entanglement pumping with single 
selection [isl Il9|. where one noisy MES is used as the 
ancilla. Here the target pair is untouched, (iii) The tar- 
get pair is purified by reducing its Z error via level-2 
entanglement pumping with double selection, where the 
successful output state of (ii) and one noisy MES are used 
as the ancillae. Here the increase in the X error is kept 
small thanks to the level-1 pumping in step (ii). At each 
step, we repeat the pumping process from a few to several 
times. Then, we finally obtain the level-2 double-pumped 
MES of high fidelity. 

The noisy MES pin, which is shared by the noisy quan- 
tum channel, is described as 



^ Fi(cr, ® cro)|0)(0|(cr, ® (To) 
i=0 



(1) 



where |0) = (|00) + |01) + |10) - |ll))/2. The pu- 
rification processes with single and double selections in 
the above pumping scheme can be described as maps 
M and F{ = Dl^'F^^'^Fl'^Fl'^/pd, re- 



Fl = S^^F^F^ 



spectively. Here, F/ and F-^"^ are the output and the 
nth input fidelities, respectively, and Ps and pd are the 
normalization factors (i.e. F^ = 1) meaning the suc- 
cess probabilities of the purification processes. Starting 
from the channel fidelity (F = Fq, Fi, F2, F3), the fidelity 
of the level-2 double-pumped MES (F = Fo,Fi,^2,^3) 
is calculated by using the tensors and D\^^ ^ which 
can be written in terms of pij and pm (see Appendix A) . 
The contour plots of the output infidelities 1 — F with 
respect to the channel fidelity F and the error probabil- 
ity pg of the local operations are shown in Fig. [3| (a), 
where Fi,2,3 = (1 - F)/?,, pij = Pg/15 and pm = Pg 
are adopted specifically. Compared to the single pump- 
ing fldi. |25||. the output fidelity is significantly improved 
in the present scheme with double selection, without in- 



creasing the number of spatial resources. 

By using the purified MES of high fidelity, we perform 
the teleportation-based two-qubit gate (TTG) between 
the data qubits stored in the nearest-neighbour QALUs. 
Including the memory space for the data qubit, a total of 
four qubits (i.e. three auxiliary qubits for the entangle- 
ment pumping and one data qubit) are required in each 
QALU. (The four-qubit system seems to be the mini- 
mum for our purpose, since we cannot purify both bit- 
and phase-flip noise by using only two auxiliary qubits.) 
The probabilities pij of the ai (g) aj errors after the TTG 
can be calculated in terms of F^, pg {pij = pg/15), and 
Pm- Roughly speaking, the gate infidelity is given by 
J2{ij)^io,o)Pij - (1 - ^) + 2pg + 2pM, which can be 
understood from the fact that the TTG is implemented 
with one purified MES, two local two-qubit gates and 
measurements [2l| (see Appendix B for details). 



IV. TOPOLOGICAL FAULT-TOLERANT 
QUANTUM COMPUTATION 

As shown above, we can perform two-qubit gates of 
high accuracy between the nearest-neighbour two QALUs 
in the teleportation-based way, but they are still sub- 
ject to small error. In order to handle this and achieve 
gate operations of arbitrary accuracy, we perform fault- 
tolerant TQC 0, where the surface code protects quan- 
tum information by virtue of the topological degeneracy 
Hi]. The logical CNOT gate operations in TQC are im- 
plemented by braiding anyons (defects) on the surface. 
Universal fault-tolerant quantum computation can be 
achieved by using the state injection and the magic state 
distillation [23 . Actually, these operations for fault- 
tolerant universal quantum computation can be executed 
by using only 2D nearest-neighbour two-qubit gates and 
single-qubit measurements Q, and therefore TQC is best 
suited for the present distributed architecture. 

The noise threshold of TQC is determined from the 
error correction procedures, which are implemented re- 
peatedly in the bulk region of the surface code ji^ . The 
syndrome measurement for the topological error correc- 
tion can be executed as shown in Fig. [3| (a) and (b), 
where the QALUs of gray and blue (red) squares indi- 
cate the data and X (Z) syndrome QALUs, respectively. 
After performing the TTGs, the blue qubits are mea- 
sured in the X basis, whose outcome corresponds to the 
eigenvalue of the stabilizer operator X'^'^ and is used to 
correct Z errors [see Fig. [11(b)]. The X errors are trans- 
formed to Z errors by the Hadamard operations, which 
is incorporated into the TTGs, and corrected in the same 
way at the next step by using the Z syndrome qubits (red 
boxes) . 

The errors during the error correction can be char- 
acterised by the independent and correlated errors on 
the three-dimensional (i.e. space-like 2D and time-like 
ID) lattice JS\. In Ref. [6], the threshold values of the 
independent and correlated error probabilities gind and 
gcor, respectively, have been numerically estimated as 
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FIG. 3: Infidelities of the purified MES and threshold curve, (a) The contour plot of the output infidelities 1 — F with respect to 
the channel fidelity F and the error probability pg of local operations, where pij — Pg/15 and pM — pg are specifically adopted, 
(b) The threshold curve with respect to the channel fidelity F and the error probability of the local operations pg — pm. 



((Zindi^cor) = (2.3%, 0.40%), whlch are attributed to the 
preparation, gate, measurement errors with equal prob- 
ability 0.75%. In om case, these probabilities gind and 
gcor a-re calculated in terms of -^2,31 Pg and pm as (see 
Appendix C for the details) 

40 

gind = i{F2+F3) + —pg+pM, (2) 
^cor = Y^Pg+PM- (3) 

Thus if they satisfy gind < 2.3% and gcor < 0.40%, fault- 
tolerance of the present architecture is ensured (the true 
threshold will be slightly higher than the value calculated 
from the above condition, which can be determined by a 
full numerical simulation) . By using these conditions and 
the output fidelity fi,2,3 obtained in the previous section, 
we calculate the fault-tolerant region of the channel fi- 
delity F and the error probability Pg of local operations, 
where pm = Pg is taken for simplicity. The threshold 
curve for these physical parameters {pg,F) is plotted in 
Fig. [3] (b). Specifically, with F ~ 1, the threshold value 
of local operations is obtained as Pr; = 0.26% {pg — 0.5% 
when Pm — 4:Pg/15 is adopted [3, 0]). It is also seen 
that the present architecture works even with very noisy 
quantum channels of fidelity F ^ 0.7 provided the error 
probability of the local operations are reasonably small 
Pg ^ 0.1%. This result contrasts with the requirements, 
F ^ 0.95 and Pg ~ lO"'^, in the ID architecture with 
five-qubit distributed systems 

The computational overhead required in the present 
architecture can be quantified by the total amount R of 
local operations plus quantum communication (i.e. the 
number of noisy MESs). The overhead R is given by 
R = KT, where T indicates the number of two-qubit 
gates in TQC, and K is the total amount of local op- 
erations plus quantum communication (i.e. the number 
of initial MESs) per TTG. Specifically, the overhead K 
is a few tens (a few hundreds) when F 0.9 (0.7) and 
Pg 0.1%, which mainly depends on the channel fidelity 



F. In order to factorise an n-bit composite number, we 
need 40n'^ Toffoli gates [l^ , each of which is implemented 
by using seven tt/S gates. Thus a total of 300n'^ 7r/8 gates 
are required to factorise a rt-bit composite number. This 
means that each tt/8 gate has to work with an error prob- 
ability of - l/(300n3). When n = 1024, such a 7r/8 gate 
requires 2 x lO^*' physical two-qubit gates under the phys- 
ical error probability at 1 /3 of the topological threshold, 
and therefore a total of T ~ 6 x 10^^ two-qubit gates are 
required in TQC. On the other hand, if 40 is required 
to achieve 1/3 of the topological threshold with F ~ 0.9 
and Pg ~ 0.1% (see Appendix D). As a result, the to- 
tal overhead for the present architecture amounts to be 
R^KT^2x 10^3. 

Although the above quantum overhead is large, com- 
parable complexity is achieved in today's classical com- 
puter within 10^ sec (a few tens of hours), where 7 x 10^ 
transistors integrated in a central processing unit (CPU) 
work at rates of 3 x lO^Hz [28] (i.e. 2 x 10^^ ~ 7 x 10* x 
3 x 10^ X 10^). It contrasts with the fact that the recent 
factorization of a general 768-bit composite number has 
taken an overhead equivalent to 1677 years with a CPU 
of 2.2GHz [H. 



V. PHYSICAL IMPLEMENTATION 

All key technologies required in the present architec- 
ture have already been demonstrated experimentally in 
various physical systems, such as trapped ions Q and 
NV-centres in diamond 0. Among them, state-of-the- 
art technologies in the trapped ion systems have already 
achieved highly accurate controls. Multiple two-qubit 
gates, single qubit rotations, and readouts have been 
achieved with error probabilities 7 x 10"^ [30l|, 2 x 10~^ 
and lO"'^ [H, [111, respectively. These benchmarks 
are comparable to the requirements for the local opera- 
tions in the QALUs. In order to attain scalability, the 
ion qubits could be distributed over separate trap zones. 
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FIG. 4: Syndrome measurements for topological error correction, (a) The syndrome measurement for topological error 
correction. The TTGs are impelented in numerical order (1— >-2— >3^4). The single line which connects the QALUs 
indicates the teleportation-based CZ gate. The double line indicates the teleportation-based CZ gate preceded by the Hadamard 
operation, where the Hadamard operation is incorporated in the TTG (see suplemental material). By measuring the qubit on 
the blue QALU, one can obtain the eigenvalue of the stabilizer operator X®*. (b) The circuit diagram in 3D (space-like 2D 
and time-like ID) of the syndrome measurement. 



There are two main schemes, which can be used as the 
quantum channels, to entangle separately trapped ion 
qubits. One is based on the quantum charge-coupled de- 
vice 34] with the microfabricated ion trap technologies, 
where entangling operations are implemented by shut- 
tling ion qubits between the storage and interaction re- 
gions. The transport of ion qubits has succeeded exper- 
imentally with high accuracy, and a two-qubit gate has 
been demonstrated with high fidelity ^ 0.9 [H, [s^. An- 
other approach utilises photons as flying qubits in order 
to entangle two separate ion qubits, where the fidelity 
~ 0.9 has been achieved experimentally jSTj . In order 
to couple ion qubits efiiciently with photons, surface- 
electrode ion traps integrated with microscale optics have 
also been investigated recently [H, [s^ . These fidelities 
of the above two entangling operations are well above our 
requirement, [pg — 5 x 10^^, pM — 10~^ and F = 0.9 
leads to the values of q^nd and gcor below the thresholds.) 

Solid state systems are another promising candidates 
for a physical implementation of the present architecture. 
In the NV-centre diamond systems, particularly, there 
are well-defined local systems as the QALUs [1, |40| • The 
control of the 12-16 dimensional systems (electric spin-1 
and two or three nuclear spin-1/2) 0, |4l| and the single- 
shot readout [i^ have already been achieved. The 
entanglement between the polarization of a single optical 
photon and a single electronic spin in a NV-centre dia- 
mond ^44] . and two-photon interference from separated 
NV-centers [H, |4^ have also been observed recently. 
These technology can be used as the quantum channel 
for the present distributed architecture. Furthermore, 



universal dynamical decoupling [47| could be employed 
to suppress noise and achieve reasonably accurate local 
operations required for the present architecture. 



VI. DISCUSSION AND CONCLUSION 

We have proposed a distributed architecture for scal- 
able quantum computation. The present architecture 
works with reasonably accurate small (four-qubit) local 
systems, which are connected by very noisy quantum 
channels. All key ingredients employed in the present 
architecture have already been experimentally demon- 
strated, and the accuracies required for them are com- 
parable to the recent experimental achievements. Ac- 
tually this is the first proposal of a practically scalable 
architecture which works well even with the preexisting 
level of the experimental devices. We believe that this 
proposal fills the gap between top-down and bottom-up 
approaches towards practically scalable quantum compu- 
tation and gives a good guideline and benchmark in the 
development of quantum devices. 

Finally let us mention that there are a lot of rooms 
to improve the performance of the present architecture. 
Instead of the bipartite entangled states, we can also 
use multipartite entangled states as resource states for 
teleportation-based multi-qubit gates. Since a part of er- 
rors during local operations, which are to be performed in 
the future computation, are removed beforehand through 
purification, the physical threshold would be improved. 
Furthermore, by using improved decoding algorithm, we 
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can fully utilise the potential power of the surface code. It 
boosts the topological threshold from 0.75% to 1.1-1.4% 
[48l |. The topological colour codes 49], whose threshold 
value 4.8% 1501 i s higher than that (3.2%) of the Kitaev's 
surface code [2g] , can be also implemented on the present 
architecture. These upgrades will improve the threshold 
values of the local operations to ^ 1%, which will become 
comparable to that of the non-distributed fault-tolerant 
schemes without modular scalability 0-01 • 
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Appendix A: Entanglement pumping 
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FIG. 5: (a) Entanglement purification with single selection, 
(b) Entanglement purification with double selection. 



the level- 1 single-pumped MES as the ancilla, where the 
phase-flip errors are checked. The purification map 1S2 
can be described as 



1. Entanglement pumping with single selection 

In the entanglement purification with single selection, 
one noisy MES is purified by using one noisy MES as the 
ancilla as shown in Fig. 2 (a). Assuming that all errors 
are given by probabilistic Pauli errors, the purification 
map can be described as 



(Al) 



where F^. indicate the fidelities of the output state. The 
transition probability tensor S]^ is given by 



(=0,3 



(A2) 



where M° and G^^^ correspond to the bilateral measure- 
ment error and bilateral CNOT gate followed by the two- 
qubit gate errors [20| . respectively, and the summation 
Si=o 3 means the postselection according to the bilateral 
measurement outcomes as shown in Fig. [5] (a), which are 
designed to check bit-flip errors. Throughout the level- 
1 single pumping, the initial noisy MES is used as the 
ancilla for pumping, and hence we take F 



(2) 



g , , ) ' Theu, the puriflcation can be 
viewed as a map from R"* ^> M*, F' = 5i(F(i)). Starting 
with the initial noisy MES (i.e. F^^' = F'"'), the single 
pumping iSi is repeatedly applied ni times. As a result, 
we obtain the level-1 single-pumped MES of the fldelity 



;iLvl 



Lvl 



(A3) 



where plvI = Yl\=o[^i^ {^™^)]k indicates the net success 
probability of the level-1 single pumping. 

In the level-2 single pumping [l^, the output state of 
the level-1 single pumping is further purified by using 



(A4) 



wher 



ere S'^ = H^SfHlH^, with i7j being the bilateral 
Hadamard transformation. By repeatedly applying the 
level-2 single pumping 712 times, we obtain the level-2 
single-pumped MES of the fidelity 



J^k^c J^aJ^b' 



7\Lv2 cn2 



= 52"^(FLvl)/pLv2 



(A5) 



The contour of the infidelity 1 — F = 10^^ of the 
level-2 single-pumped MES is plotted with respect 
to the channel fidelity F and error probability 
Pg — Pm of local operations for each (ni,n2) = 
(2, 4), (3, 4), (3, 7), (5, 6), (5, 8), (5, 10), (5, 11), (5, 13) in 
Fig. [5] (a). There is a tradeoff between the channel 
fidelity F and error probability Pg = pm of local 
operations when the the numbers of repetitions (ni, 712) 
are increased. This behavior can be understood that 
either bit or phase error is checked repeatedly at each 
pumping level, which is designed to obtain the output 
MES of high fidelity even with low channel fidelity F as 
discussed in Ref. [16]. 



2. Entanglement pumping with double selection 

In the level-1 entanglement pumping with double se- 
lection, one noisy MES is purified by using two initial 
noisy MESs as the ancillae as shown in Fig. [S] (b). Sim- 
ilarly to the previous case, the purification map Pi can 
be described as 



E 



T-)ijk 771(1) ciini T^ini 
Iran i j k i 



(A6) 



m=0,3;n=0,l 



where the transition probability tensor is given by 
Dl^n = HtM^M^iC^GZ with = HiM^Hf. Start- 
ing with the initial noisy MES (i.e. F^^) = F™'), the 
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The contour of the infidehty 1 — F = 10~^ is plotted with 
respect to the channel fidelity F and error probability 
Pg = Pm of the local operations for each (ni, rrii, 7712) = 
(2, 5, 5), (2, 4, 8), (3, 3, 9), (3, 3, 11), (3, 3, 13), (3, 4, 14) in 
Fig. |6](b). In Fig. 3 (a), the result with (711,7711,7712) = 
(3, 4, 14) is plotted. In the level-2 double pumping, the 
acceptable range of the channel fidelity F is improved 
without sacrificing that of the error probability pg of local 
operations, when the numbers of repetitions (ni, 1711,1712) 
are increased. This is due to the following two factors: (i) 
Both bit and phase errors are checked at each pumping 
level, (ii) A large amount of operational errors, which 
are left on the output state of single selection, can be 
detected by double selection [20| . 
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FIG. 6: Infidelities of the purified MESs. (a) The contour of 
the infidelity l-F = 10"^ of the level-2 single-pumped MES 
is plotted with respect to the channel fidelity F and error 
probability pg = pM of local operations for each (711,712) = 
(2, 4), (3, 4), (3, 7), (5, 6), (5, 8), (5, 10),_(5, 11), (5, 13). (b) The 
contour of the infidelity 1 — F = 10 ^ of the 
level-2 double-pumped MES is plotted with respect to 
the channel fidelity F and error probability pg = 
Pm of the local operations for each (711,7711,7772) ~ 
(2, 5, 5), (2, 4, 8), (3, 3, 9), (3, 3, 11), (3, 3, 13), (3, 4, 14). 



level-1 double pumping Vi is repeatedly applied tui 
times, and then we obtain the level-1 double-pumped 
MES of the fidelity 



Pr^(F"")/rLvi 



(A7) 



where tlvI is the net success probability similarly to the 
previous case with single selection. In the level-2 double 
pumping, the output state of the level-1 double pump- 
ing is further purified by using the level-1 single pumped 
MES and initial noisy MES as the ancillae. The level-2 
double pumping ^2 can be described as 



F' = 



m=0,3;ri=0,l 



TTiini 



(A8) 



By applying 2?2 repeatedly m2 times, we obtain the level- 
2 double-pumped MES of the fidelity 



pLv2 =p™2(pLvl)/^^ 



v2 



(A9) 



Appendix B: Error analysis of TTGs 

We employ three types of TTGs in the syndrome mea- 
surement (see Fig. 4 (a)). 



1. TTG of type I 

The TTG of type I [® in Fig. 4 (a)] is as follows: 



syndrome qubit 

1+) 



data 
qubit 



H 



purified MES 



— •- 



—j ^x^" 



(Bl) 



where the Hadamard operations is required to transform 
Z errors to X errors. The probabilities of ai (g) aj errors 
on the two output qubits are given by 



Ai 


I 


X 


Y 


z 


I 




PX{I,X} 


Py{i,x} 


Fi +Pz{i,x} 


X 










Y 










Z 


F3 +PX{Y,Z} 


PM +Pl{Y,Z} 


PZ{Y,Z} 


F2 + Py{y,z} 



where we switched the notation from pij to pab {A, B G 
{I,X,Y,Z}), and Pa{b,c} means pab +Pac- 



2. TTG of type II 

The TTG of type II [® and ® in Fig. 4 (a)] is as 
follows: 



syndrome qubit 



e — ^2 



data 
qubit 




(B2) 
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The probabilities of Ui ® (Jj errors on the two output 
qubits are given by 



data qubit 



Ai 


I 


X 


Y 


z 


I 




PX{I,Z} 


Py{i.z} 


Fi +pz{i,z} 

+P\x,Y}X 


X 


P\i,z\x 






PM + P'{X,Y}I 


Y 


P{I,Z}Y 






P{X,Y}Z 


Z 


F3 +PX{X.Y} 
+P'{I,Z}Z 


Pm +Pi{x,y} 


PZ{X,Y} 


F2 +Py{x,y} 



where pij and p'^j indicate the error probabihties of the 
top and bottom two-qubit gates, respectively, as shown 
in the above diagram. 



3. TTG of type III 

The TTG of type III [(3) in Fig. 4 (a)] is as follows: 



syndrome qubit 



data 
qubit 



H 



purified MES 



{Fi, F2, F3} - 
— •- ' 



{Vii} 



{Pij 



. Pm 



(B3) 



The probabilities of Ui ® aj errors on the two output 
qubits are given by 





I 


X 


Y 


z 


I 




Px{i,x} 


Py{i,x} 


Fi +Pz(i,x} 

'^P\x,Y}X 


X 


P'{I,Z}X 






PM + P\x,Y}I 


Y 


P'{I,Z}Y 






P{X,Y}Z 


Z 


F:i + Px{Y,z} 

+P'{I,Z}Z 


PM +Pl{Y,Z} 


PZ{Y,Z} 


F2 + Py{y,z} 

+P'{X,Y}Y 



where pij and p[j indicate the error probabilities of the 
top (CNOT) and bottom (CZ) two-qubit gates, respec- 
tively, as shown in the above diagram. 



Appendix C: Threshold analysis 

The syndrome measurements for topological error cor- 
rection are implemented as shown in Fig. [T] The errors 
during the error correction on the primal lattice (6l, [5lll52j 
can be described as Z errors located at a, b, and c, which 
are depicted with green lines in Fig. [T] The errors at a 
and 6, c correspond to errors on the measured syndrome 
and data qubits, respectively. Similarly, the errors dur- 
ing the error correction on the dual lattice are located 
at a, 5, and c, which are depicted with orange lines in 
Fig. [71 Since the error corrections on the primal and 
dual lattices are implemented independently, we can ig- 
nore the correlations between errors on the primal and 




syndrome qubit 



FIG. 7: The unit cell of the syndrome measurements for topo- 
logical error correction. The two-qubit gates are implemented 
in the numerical order. The errors during the error correction 
on the primal and dual lattices can be described as Z errors 
located at a, b, c (green lines) and a, b, c (orange lines), respec- 
tively. 



dual lattices (i.e. correlated errors between a,6, c and 
a,b,c). Furthermore, behavior of the errors on the dual 
lattice is equivalent to that on the primal lattice, since 
they are symmetric under reflection and rotation. Thus 
we consider only the errors on the primal lattice, i.e., 
those errors at a, b, and c. 

Let g", q'', and q'^ denote the probabilities of the in- 
dependent Z errors at a, 6, and c, respectively. If er- 
rors are located on the data qubits independently with 
equal probability, and syndrome measurements are per- 
fect, then g" = 0, q'^ = q"^. In such a case, the threshold 
value is given by q'' = q" = 11% With q" = q^ ^ q'= 
(i.e. independent errors on the syndrome and data qubits 
with equal probability), the threshold values have been 
estimated as q'' = q'> = q'= = 2.9% ^ and 3.3% [H by 
using the minimum-weight perfect matching algorithm 
and the random-plaquette Z2 gauge theory, respectively. 

In the present case (and also in Ref. [1, HI]): however, 
we have to take the correlated errors into account, since 
the two-qubit gates for the syndrome measurements in- 
troduce correlation. By using the commutation relations 
between errors and two-qubit gates, it is found that the 
correlated errors are located only at (a, b) and (a, c) of 
the same unit cell and {b,b) of the neighboring unit cells. 
Such probabilities are denoted by g"-'', g"^^, and q''^''. The 
probabilities of independent and correlated errors can be 
calculated in terms of the error probabilities p^'^ of the 
Ith two-qubit gate (see Fig. [7]), where A and B indicate 
errors on the syndrome and data qubits respectively, and 
the Hadamard operation is included in the two-qubit gate 
for I = 1,3,5,7: 



= + + £J + Pi'i +PP+ PM. (CI) 
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^P'jj+p'il (C2) 








^pW+I^, (C3) 




(7) 
= Pzx - 


ypfl 


(C4) 


qa.c 


(5) 
= Pzx - 




(C5) 




(2) 
= Pxi - 


^pI^- 


(C6) 



tions (IC7I1 read 



where pp +pm in 9° corresponds to the preparation and 
measurement errors of the ancilla qubit for the syndrome 
measurement, and the subscripts z, z , x, and x mean 
J2a=y,z^ Y.a=i.x^ Y.a=x.y^ and Y.a=i..z respectively. 

For example, p^'j = Y.a=y,zT.b=i,zP% ^ Pyi+Pz] + 

(0 I (0 
Pyz + Pzz ■ 

In the original scheme [6] , two-qubit gate errors B 
{A, B = I,X, Y, Z) and single-qubit gate errors A {A ^ 
X, Y, Z) are assumed to occur with equal probability 
Pg/15 and Pg/3, respectively. Thus p^'^ = Pg/15 for all 



A, B when I = 2, A, 6, 8, and p\% = = p'zy = 6pg/15 

and = Pff/15 for other A,B when Z = 1,3,5,7. 

The preparation and measurement error probabilities are 
taken as pp = pm = Pg- Then, Eqs. (|Cip - (IC6p read 
q'' = 46pg/15, q'' = Upg/15, q" = Upg/15, q"''' = q""" = 
q^'^ — 8pg/15. On the other hand^the threshold value for 
Pg is given by Pg = 0.75% in Ref . Q , which is obtained by 
using the minimum-weight-perfect-matching algorithm. 
This leads to threshold conditions of these probabilities 
(sufficient conditions for fault-tolerance) as 



q'' < 0.023, q'' ^q" < 0.022, 



q^-^ < 0.0040. 



(C7) 



In the present architecture, on the other hand, we use 
the TTGs of type I (for I = 1,5), II (for / = 2,4,6,8), 
and III (for I = 3,7), where the state preparation for 
the syndrome measurement and Hadamard operations 
are incorporated (see TTG I and III). Specifically, for 
the TTG of type I (i.e. I = 1,5), the error probabilities 
Pal ('^j b = z,z, X, x) are given by 



= 4pg/15+pM, 
pi2 = F2 + ^3 + 4pg/15, 

Pzx 

H) 



4pg/15, 



ptl = Fi F2 8pg/15, 



(C8) 
(C9) 
(CIO) 
(Cll) 
(C12) 



where we take pab = Pgl^^- Similarly, for the TTG of 
type II and III (i.e. I ^ 1,5), they are given by 

= p'il^Apg/lb+PM, 



(C13) 
(C14) 
(C15) 

p'ii - Fi + + 12pg/15. (C16) 
By using the these, Eqs. (|Cip - (|C6p and threshold condi- 



Pzx 

P^l = F2+h + \2pg/\b, 

Psl = Pxl = 4pg/15, 



40 

g"^ = 4(F2 + Fs) + —Pg + PM < 0.023, (C17) 

q'' - 2(A+A) + ^Pg + 2pM < 0.022, (C18) 
1 o 

32 

q^ = 2(Fi + Fa) + — Pg + 2pM < 0.022, (C19) 

1 D 



g"'" = (z'''' = —Pg+PM< 0.0040. (C20) 



(Note that pp = in the present case, since prepa- 
ration of |-|-) ancilla for the syndrome measurement is 
taken into account in the TTG of type I.) The thresh- 
old curve {F,pg = pm) are plotted in Fig. 3 (b), 
which is calculated with F obtained by using the tran- 
sition probability tensors. Specifically, in the limit of 
F — > 1 , the output fidelity of the purification is obtained 
as (Fi^2,F3) = (4pg/15,2pg/15,2pg/15) in the leading 
order [20]. With pM — Pg, we obtain the following con- 
ditions on Pg-. 



71 

g° = —p„ < 0.023 < 0.0049, 
15 

82 

q'' ^ —pn < 0.022 4^ p„ < 0.0040, 
15 



74 
15 



Pg < 0.022 ^ pg < 0.0045, 



(C21) 
(C22) 
(C23) 



23 



—Pn < 0.0040 ^Pn < 0.0026. 
15 

(C24) 



This leads to the physical threshold Pg = 0.26% with 
F ~ 1. Since q"-''^^'^ are smaller than the threshold values 
for them, the true physical threshold would be higher 
than Pg = 0.26%, which would be determined by using a 
full numerical simulation. 



Appendix D: Resource analysis 

The operational overhead under the physical error 
probability at 1/3 of the topological threshold is plot- 
ted as a function of the circuit size in Fig. 11 of Ref. 
It reads that if we want to perform J7 = 3 x 10^^ 7r/8 gates 
accurately (i.e. an accuracy of ^ is required for each 
7r/8 gate), a logical 7r/8 gate requires T = 2 x 10^° (O3 
in Ref. [6|) physical two-qubit gates. 

In Fig. [8] (a), 1/3 of the topological thresholds in 
the present case [i.e., the curves {F,pg) which satisfy 
q'' < 0.0049/3, q'' < 0.0040/3, q'= < 0.0045/3 and 
ga,b ^ ^ qbM < 0.040/3] are plotted for each 

(ni, mi, 7712) - (1,1, 2), (1,2,2), (1,2,3), (1,2,4). It can 
be seen that {ni, 1711,1112) = (1,2,2) achieves 1/3 of the 
topological threshold even with F ~ 0.9 and pg ^ 0.1%. 
The the total amount K of local operations plus quan- 
tum communication (i.e. the number of initial MESs) 
per TTG (i.e. per purified MES) is determined for given 
repetition numbers (ni, 7711,7772) by the channel fidelity F 
and error probability Pg = pm of local operations through 
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FIG. 8: {F,pg) at 1/3 of the topological threshold and operational overheads per TTG. (a) The channel fidelity F and error 
probability pg = pm of local operations with which the independent and correlated error probabilities are at most 1 /3 of the 
topological threshold are plotted for each (ni,mi,m2) = (1, 1, 2), (1, 2, 1), (1, 2, 2), (1, 2, 3), (1, 2, 4). (b) The contours of the 
total amount K of local operations plus quantum communication (i.e. the number of initial MESs) per TTG are plotted with 
respect to the channel fidelity F and error probability pg — pm of local operations for (ni, mi, 7712) — (1, 2, 2). 



the success probabilities plvI, ''LvI, and rLv2- In Fig. [5] Pg = Pm- Specifically, if ~ 40 when F ^ 0.9 and Pg ~ 
(b), the contours K = 30,60,120 of the total amount 0.1%. Accordingly, the total operational overhead can be 
K for (ni,mi,m2) = (1,2,2) are plotted against F and calculated as i? = KT - 40x2xl0i°x3xl0" - 2x10^3. 
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